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Abstract. We considered the structure of steady-state plane-parallel radiative shock waves propagating through the partially 
ionized hydrogen gas of temperature T\ = 3000 K and density 10~ 12 gm cirT 3 < p\ < 10~ 9 gm cm 4 . The upstream Mach 
numbers range within 6 < M\ < 14. In frequency intervals of hydrogen lines the radiation field was treated using the transfer 
equation in the frame of the observer for the moving medium, whereas the continuum radiation was calculated for the static 
medium. Doppler shifts in Balmer emission lines of the radiation flux emerging from the upstream boundary of the shock wave 
model were found to be roughly one-third of the shock wave velocity: —6V « \U\. The gas emitting the Balmer line radiation 
is located at the rear of the shock wave in the hydrogen recombination zone where the gas flow velocity in the frame of the 
observer is approximately one-half of the shock wave velocity: -V* ~ hU\. The ratio of the Doppler shift to the gas flow 
velocity of SV/V* ~ 0.7 results both from the small optical thickness of the shock wave in line frequencies and the anisotropy 
of the radiation field typical for the slab geometry. In the ambient gas with density of p\ ^ 10~" gm cirT 3 the flux in the Ho- 
frequency interval reveals the double structure of the profile. A weaker H/? profile doubling was found for p, > 10~ 10 gm cirr 3 
and U\ < 50 km s . The unshifted redward component of the double profile is due to photodeexcitation accompanying the 
rapid growth of collisional ionization in the narrow layer in front of the discontinuous jump. 
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1. Introduction 

It is now a well-established fact that hydrogen emission lines 
are a characteristic feature of radially pulsating stars of vari- 
ous t ypes. Strong hydrogen emission is observed in Mira type 



ous types. 5>t 
J.Tovlll947[[i 



1954), W Vir (Abt 1954; Wallerstein 1959}) and 



RV Tau dPreston 1 1962 ) pulsating variables. Moreover, in the 
spectra of Mira stars Balmer emis sion lines persist during the 
major part of the pulsation period (iJovll 19471 iRichter & W ood 
2001). The intensity of the hydrogen emission seems to corre- 
late with the amplitude of the puls ation sin ce RR Lyr variables 
exhibit only weak emission lines (Preston 1964) and in classi- 
cal Cepheids (for example i n Dor) hydrogen emission lines 
are scarcely detected ( Hutc hinson et alj|l975h . The hydrogen 
emission in the spectra of pulsating stars is thought to be due 
to radiative cooling of the gas compressed by the shock wave 
propagating throu gh the stellar at mospher e in each pulsation 
cvcle JKraft et alll959tlWailersteirJl959tlAbt & Hardiell960t 
lr,orbatskiill96lh . 

High resolution spectroscopy reveals the doubling of Ha 
and H/? emission profiles, whereas profiles of higher Balmer 
lines exhibit only the asymmetry. This feature is observed not 
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only in Mira stars jBidelman & Ratcliff el 19541 Ifox et alJl984l 
lGilletlll988HWoodsworthl 19951) but als o in W Vir and RV Tau 
variab les JLebre & Gilleflll99ll Il992l) . iBidelman & Ratcliffd 
il 19541) explained the Ha profile doubling observed in the Mira 
type star T Cen as an absorption reversal rather than a real 
duplicity, the absorption resulting from the cool hy drogen gas 
above the propagating shock wave. IWillsonl J 19761) interpreted 
the double structure of the emission profiles in terms of a spher- 
ically symmetric shock wave with a radial distance from the 
center of the star appreciably larger than the radius of the pho- 
tosphere. According to this model the flux of the redshifted 
component emerges from the bac k side of the shock wave 
moving outward from the observer. Woo dsworthl i 19951) mod- 
elled double Ha profiles as three emission components of equal 
width, two of which are blended. Thus, different phenomeno- 
logical models demonstrate the ambiguity existing so far in our 
understanding of the origin of the Balmer emission lines pro- 
duced by the shock waves in stellar atmospheres. 

This paper is the fifth in our series on the struc- 
ture of radiative shock waves. In our previous Papers I -IV 
llFadevev & Gilletl fl998. 20001 1200 it iFadevev et al]l2002l) we 
presented the method of computation and described results for 
the structure of radiative shock waves propagating through par- 
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tially ionized hydrogen gas with temperature and density typi- 
cal for atmospheres of pulsating late-type stars. An advantage 
of this approach is that the gas dynamics, radiation field and 
atomic level populations are considered self-consistently for 
the whole shock wave model. However, in all our previous pa- 
pers the radiative transfer was treated in a static medium ap- 
proximation, so that we were unable to compare the calculated 
monochromatic radiation flux with observed emission profiles. 
Below we describe the shock wave models with Doppler shifts 
in the line profiles computed from the transfer equation in the 
frame of the observer. 

2. Method of computation 

As in our previous studies we consider a steady-state plane- 
parallel shock wave propagating through a homogeneous 
medium consisting solely of atomic hydrogen gas. The shock 
wave model is represented by a flat finite slab moving together 
with a discontinuous jump toward the observer at velocity V s 
which is regarded as negative. The unperturbed gas enters the 
upstream face of the slab with velocity U\ = — V s , and through- 
out this paper U\ is referred as the shock wave velocity. 

The space coordinate X is measured along the path of the 
gas element moving through the slab and is within the range 
X\ < X < Xn, where X\ and X^ are the coordinates of the up- 
stream and downstream faces of the slab, respectively. We set 
X = at the discontinuous jump where the gas flow velocity 
U and the gas density p undergo an abrupt change. At the up- 
stream face of the slab the gas is assumed to be unperturbed, 
whereas at the downstream slab boundary the electron temper- 
ature T e and the hydrogen ionization degree xh are assumed 
to approach their postshock asymptotic values. In general, the 
space coordinates of both boundaries were determined from 
trial calculations in such a way that the solution describing the 
structure of the shock wave is independent of coordinates X\ 
and X?j. 

For the calcula tion of the shock wave structure we follow 
ourprevios studies (Paper I; Paper II) and solve the initial value 
problem for the fluid dynamics and rate equations treated as a 
system of ordinary differential equations, while radiative trans- 
fer is treated as a two-point boundary value problem. To take 
into account the strong coupling between the gas flow and the 
radiation field produced by this flow the ordinary differential 
equations and the transfer equation are solved iteratively. 

The radiation transfer equatio n is solved with the improved 
Feautrier method ( Rvbicki & Hummer 1991) and the slab is di- 
vided into N = 1.6 x 10 4 cells with 6000 cells ahead of the 
discontinuous jump. Space intervals AXj-\/2 = Xj - Xj-y in- 
crease in a geometrical progression in both directions from the 
cell interface Xg = 0. For the sake of convenience the vari- 
ables at the cell centers nearest to the discontinuous jump are 
denoted with superscripts minus and plus, that is X~ = Xg_i/2 
andX + = X a+ i/2. 

The frequency range 7 x 10 13 Hz < v < 10 16 Hz consists of 
a set of intervals with boundaries either at hydrogen ionization 
thresholds or at limits of line frequency intervals. In the present 
study we consider the hydrogen atom with L = 5 bound levels 
and the continuum, so that the full frequency range is repre- 



sented by 15 intervals for the continuum and by 10 intervals 
for lines. 

The electron collision excitati on rates were c alcu lated with 
effective collision str engths from ISchplz et all Jl990l) for the 
1 — > 2 transition, from Callawav ( 1994) for the 1 — > 3 and 2 — > 
3 transitions and from Aggarwal et al. (1991) for other bound- 
bound transitions. The electron coll ision ionization rates were 
computed with polynomial fits from lLennon et alJ fl9 86i The 
photoexcitation rates were calculated for Doppler profiles. 

In general, the interaction of the radiation field with the 
moving medium is described in terms of the comoving frame 
partial differential transfer equation. The method developed by 
iMihalas et all ill 9751) treats the comoving frame transfer equa- 
tion as a two-point boundary value problem with respect to the 
spatial coordinate and as a Cauchy problem with respect to the 
frequency. However our attempts to apply such an approach 
to the shock wave structure failed due to numerical instabil- 
ity arising because of the existence of the velocity jump. The 
difficulties accompanying the solution of the comoving frame 
transfer equation for no nmonotonic ga s flow velocity fields 
were discussed earlier by Mihalas ( 1980). 

In our study the ratio of the gas flow velocity to the speed 
of light is confined to values of U/c < 3 • 10~ 4 , so that the 
only significant effect of the moving medium is the Doppler 
shift in the spectral lines. The radiation flux within all line fre- 
quency intervals is less than ~ 10~ 3 of the total radiative flux 
produced by the shock wave, so that we assume that without 
significant loss of self-consistency the Doppler shifts can be 
regarded as insignificant for the structure of the shock wave. 
Such an assumption allowed us to calculate the shock wave 
model iteratively with a static transfer equation. When the it- 
eration procedure converges, we solve the transfer equation in 
line frequency intervals for the moving medium. In support of 
this approach below we compare solutions of the transfer equa- 
tion for a moving and for a static medium. 

In plane parallel geometry the transfer equation in the frame 
of the observer is llMihalasl 1978b 

A* ~^T- = ^vOj) -*vG")/y(M), (1) 

where ft = cos(ff) is the directional cosine, I v (ju) is the spe- 
cific intensity, rj v (fj) and Xv(M) 816 tne emission and extinction 
coefficients, respectively. For brevity in Eq. Q we omit the 
functional dependencies on X. 

The motion toward the observer is regarded as negative, 
whereas the velocity in the frame of the shock U is positive. 
Velocities in both frames relate via 

v = v- + u--u, (X< 0), 

+ + w 

V = V + + U + - U, (X> 0), 

where V~ = U~ - U\ and V + = U + - U\. Thus, the postshock 
compression ratio expressed in terms of the gas flow velocity 
in the frame of the observer is given by 

L -_ Yi . (3) 

p l V S -2V + + V 

Though on the r.h.s. of Eq. (TJ} we have only the observer's 
frame velocity, this expression is equivalent to Eq. (8) given in 
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Paper III where the postshock compression ratio is expressed in 
terms of the specific internal energy of the gas and the radiation 
flux. 

The line opacity and line emissivity coefficients k c v (h) and 
rfyifj) are calculated on the assumption of complete redistribu- 
tion, with the Doppler profile given by 



logp, = -10 T, = 3000K U, = 70 km 



0(v) = 



1 



: exp 



• - v + ^//Vj /Av 2 D 



where vo is the line-center frequency, 



v / 2kT a 

Av D = — 

c \ m H 



1/2 



(4) 



(5) 



is the Doppler width, and oth is the mass of hydrogen atom. 

3. Shock wave model 

The self-consistent solution of the equations of fluid dynam- 
ics, radiative transfer and rate equations for atomic level popu- 
lations describes the structure of the shock wave and depends 
on three general parameters, which are: the density p\ and the 
temperature T\ of the unperturbed gas, and the velocity of the 
shock wave U\ with respect to the ambient gas. The structure 
of the radiative shock wave is shown in Fig^where the follow- 
ing are plotted as a function of space coordinate X: the electron 
temperature T e and the temperature of hydrogen atoms T a (here 
we assume that neutral hydrogen atoms and hydrogen ions have 
the same temperature), the hydrogen ionization degree xe, the 
total radiative flux Fr, and the Lyman continuum flux Fj_, y . For 
the sake of convenience we use a logarithmic scale with respect 
to the space coordinate, so that each variable is represented 
by two separate plots for the preshock (X < 0) and postshock 
(X > 0) regions, respectively. 

In the partially ionized hydrogen gas the upstream bound- 
ary of the shock wave is associated with ionization of hydrogen 
atoms ahead of the discontinuous jump, whereas the down- 
stream boundary is associated with the recombination zone. 
Within the major part of the plots shown in the lower panel 
of Fig.[2the total radiative flux is upstream (i.e. Fr < 0) since 
the radiation is generated at the rear of the shock wave in the 
hydrogen recombination zone. The total radiative fluxes emerg- 
ing from both faces of the slab in opposite directions are nearly 
equal, that is, -Fri « Frn. 

The interaction of the gas flow with the radiative field due 
to absorption or emission of radiation is described in terms of 
the divergence of the total radiative flux 



V-F, 



DO 

,.4,/ 



(j]y - Kyjy) dV, 



(6) 



where J v is the mean intensity. In the upper panel of Fig. 13 we 
show the plots of V • Fr computed for the moving and for the 
static medium. In the preshock region (-10 4 cm < X < cm) 
the gas is radiatively heated due to absorption of the Lyman 
continuum radiation. Just ahead of the discontinuous jump 
(-10 2 cm < X < cm) radiative heating substantially increases 
because of absorption of radiation in Lyman lines though at 



10' 



-l 



log (-X) 



2 

(X) 



Fig. 1. The structure of the radiative shock wave. The left- and 
right-hand sides of the plots represent the preshock and post- 
shock regions, respectively. Upper panel: the electron temper- 
ature T e (solid line) and the temperature of neutral hydrogen 
atoms and hydrogen ions T a (dashed line). Middle panel: the 
hydrogen ionization degree x^. Lower panel: the total radiative 
flux Fr (solid line) and the Lyman continuum flux F]_ y (dashed 
line). The space coordinate X is in centimeters. 

Balmer frequencies the gas emits more energy than it absorbs. 
The radiative cooling of the gas due to emission in Balmer lines 
is less by an order of magnitude than radiative heating due to 
absorption in Lyman lines. 

Such a difference in behaviour of the radiation field in the 
La and Ha frequency intervals is clearly seen from the com- 
parison of the plots of the divergence of the line radiation flux 



V ■ F f = An 



vo+Av 



(ju) - k v (h)J v ] dv 



(7) 



vq-Av 



shown in the middle and lower panels of Fig. |5] Here Av = 
3.5 x 10~ 4 vo is the half-width of the line frequency interval. 
The opacity and emissivity coefficients are k v (jj.) = K e v (fj) + k° v 
and t] v (jj.) = rf v {pL) + rf v , respectively, where superscripts I and c 
denote the line and the continuum variables. Thus, expression 
defines the divergence of the radiative flux for both line and 
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logp,= -10 Tj = 3000K U, = 70 km s 




4 2 2 4 

log(-X) log(X) 



Fig. 2. The divergence of the total radiative flux (upper panel), 
of the La radiative flux (middle panel) and of the Ho- radiative 
flux (lower panel) as a function of space coordinate X. Solid 
and dashed lines show the solutions for the moving and for the 
static medium, respectively. 



continuum radiation. It should be noted, however, that the con- 
tribution of the continuum radiation can be neglected because 
the maximum values of the line opacity and line emissivity co- 
efficients dominate the continuum values by as much as two 
orders of magnitude. 

Most of the spectral line radiation is transported in Lyman 
series lines. As is seen in the middle panel of Fig. 13 Doppler 
shifts lead to stronger radiative heating of the gas ahead of 
the discontinuous jump and to stronger radiative cooling be- 
hind the discontinuous jump. However the size of this zone is 
about four orders of magnitude less than the length of the shock 
wave, so that the effects of Doppler shifts are small enough 
and can be neglected. For example, the change of the upstream 
La radiative flux due to Doppler shifts does not exceed 10%, 
whereas the La radiative flux emerging from the upstream face 
of the slab is less than 10~ 3 of the total radiative flux. For the 
model shown in Figs.^ an d|2]the optical thickness of the slab 
at La and Ha line-center frequencies is ro(La) w 6 x 10 6 and 
To(Ha) w 40, respectively. 



log p, = -10 T! = 3000K U 1 = 70 km s 




-100 100 

km s _1 



Fig. 3. The divergence of the monochromatic radiative flux 
within the La frequency interval just ahead of the discontin- 
uous jump (upper panel) and at the maximum of V ■ ¥[(La) 
(lower panel). Solid and dashed lines show the solutions for 
the moving and for the static medium, respectively. 

Thus, the effects of Doppler shifts are most pronounced in 
the close vicinity of the discontinuous jump. In Fig.[3]we show 
the plots of the divergence of the monochromatic radiative flux 
V • F v (Lq') as a function of frequency v for the cell center X~ 
just ahead the discontinuous jump and for the maximum of V • 
Ff (La) behind the discontinuous jump (see the middle panel of 
Fig. El 

4. Formation of line radiation 

Within the shock wave one can delineate three distinct zones 
with substantial changes of level populations n, due to radiative 
bound-bound transitions that lead to the appreciable emission 
of spectral line radiation. To compare the effect that different 
atomic transitions have on the line radiation in transition i — > j 
we calculate the net photoexcitation rates 

Zij^mRij-n^R] (8) 

j 

and the net collisional excitation rates 

Yjj = n e dj yn - n j(n*/n*j)\ , (9) 
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where Rjj and (n*/npl?T. are the photoexcitation and photode- 
excitation rates, respectively, Cy is the rate of collisional exci- 
tation and n*/n* is the thermodynamic equilibrium ratio given 
by the Saha-Boltzmann relation for a nonequilibrium number 
density of electrons « e . We also estimate the net photoioniza- 
tion rates 



Zik = riiRi k - n e —R],, 



(10) 



where R^ and (n*/nl)R . k are the photoionization and photore- 
combination rates, respectively, and the net collisional ioniza- 
tion rates 



Y ik = n,n e C,7t 1 - 



(11) 



where C& is the collisional ionization rate. 



4.1. The radiative precursor 

As can be seen in Fig. 0] the only essential atomic transition 
within the major part of the radiative precursor is photoioniza- 
tion from the ground state which is due to the absorption of 
the Lyman continuum radiation. However in the close vicin- 
ity of the discontinuous jump (-10 2 cm < X < 0) the rapid 
growth of collisional ionizations is accompanied by photodeex- 
citation while the rates of collisional bound-bound transitions 
remain negligible. The effect of preshock photodeexcitation on 
the structure of the shock wave is quite small; however, as we 
show below, they lead to remarkable features in the Ha and H/3 
profiles of the emergent radiation flux. 

To estimate the contribution of preshock photodeexcita- 
tions 3 — > 2 on the emergent Ha radiation flux 



v (Ha)+Av 



F;(H 



F y dv 



(12) 



VQ(Ha)-Av 



we calculated the quantity 



u 

^(Hcir) = - J 

A', 



V • F ( (Ha)dX. 



(13) 



The ratio J'iHa) / F '[{Ha) was found to decrease with decreas- 
ing gas density p\. For example, in shock waves with velocity 
U\ ~ 70 km s the fraction of the preshock Ha emission in 
the emergent radiation flux F((Ha) is within the range 10~ 2 < 
3(Ha) I F { (Ha) < 0.2 for 10~ 12 gmcirT 3 < p, < 10~ 9 gmcirT 3 . 



4.2. The zone of hydrogen ionization 

Hydrogen ionization behind the discontinuous jump proceeds 
collisionally since photoionization rates are much lower than 
the rapidly increasing collisional ionization rates (see panel (b) 
in Fig. |5}- The continuum radiation flux remains nearly con- 
stant so that any changes of the total radiation flux behind the 
discontinuous jump result mainly from La emission. In panel 
(c) in Fig. |5]we show the net collisional and photoexcitation 



0.5 



0.0 




Fig. 4. The zone of the radiative precursor. Upper panel: the 
hydrogen ionization degree. Lower panel: net collisional ion- 
ization rates (solid lines) and net photoionization and photoex- 
citation rates (dashed lines). 

rates for transitions 1-2 and 1-3. Even though photodeex- 
citation 2 — > 1 does not exactly balance collisional excitation 
1 — > 2, it is responsible for strong emission of La photons 
and leads to appreciable radiative cooling of the gas seen in the 
middle panel in Fig.|5]as the maximum of V • F((La). Bound- 
bound transitions from the ground state to levels j > 3 are 
collisionally dominated and emission of radiation due to pho- 
todeexcitation is quite small. Among excitations from levels 
i > 2 one should only note the collisional transition 2 — > 3 ac- 
companied by emission of Ha photons (see panel (d) in Fig. [5}, 
which however is scarcely perceptible. 

4.3. The zone of hydrogen recombination 

The center of the radiation emitting zone we define as a layer 
where the total radiative flux over the upstream hemisphere 
(— 1 < fJ. < 0) equals the total radiative flux over the down- 
stream hemisphere (0 < fi < 1), that is, 



oo u oo 1 

- J J Iy(ii)iid[idv = J J I y (jj)iMiiJdv. 



(14) 



o o 



Thus, in the center of the radiation emitting zone we have Fr = 
0. For the model shown in Fig. ^ condition fl!4i is fulfilled in 




Fig. 5. The zone of hydrogen ionization behind the discontin- 
uous jump, (a) The hydrogen ionization degree, (b) Net colli- 
sional ionization rates (solid lines) and the net photoionization 
rate Z\k (dashed line), (c), (d) Net collisional excitation rates 
(solid lines) and net photoexcitation rates (dashed lines). 



Fig. 6. The zone of hydrogen recombination at the rear of the 
shock wave, (a) The hydrogen ionization degree, (b) Net col- 
lisional ionization rates (solid lines) and net photoionization 
rates (dashed lines), (c), (d) Net collisional excitation rates 
(solid lines) and net photoexcitation rates (dashed lines). 



the hydrogen recombination zone at xh = 0.69. It should be 
noted, however, that for the monochromatic radiation flux F v 
the space coordinate X(F V = 0) is frequency dependent. For 
example, the radiative flux averaged over the Lyman continuum 
changes its sign closer to the discontinuous jump in layers with 
x H ~ 0.99. 

As is seen in panel (b) in Fig. [6] the initial hydrogen re- 
combination results from radiative downward transitions to the 
ground state. However, at ionization degrees jch < 0.94 recom- 



bination of hydrogen atoms is due to downward bound-free 
transitions on levels i > 2, both radiative and collisional, while 
atoms in the ground state undergo ionization. Thus, the final 
relaxation of the postshock gas is appreciably slowed down by 
ionization of the ground state hydrogen atoms. 

In the recombination zone the ground-state hydrogen 
atoms undergo collisional excitation, though radiative bound- 
bound transitions have the opposite direction. Within the re- 
combination zone the most conspicuous transition is excitation 
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1 — > 2 which is accompanied by photodeexcitation 2 — > 1 with 
resulting rate Y\2 + Zn < (see panel (c) in Fig.[6j. Though 
in the hydrogen recombination zone the absolute value of the 
net rate Z\i is about an order of magnitude less than in the ion- 
ization zone (compare plots in panel (c) in Fig. [6] with those 
in panel (c) in Fig- EJ, most La photons of the shock wave are 
created in the recombination zone because of its much larger 
length. Collisional excitation 1 — > j for j > 3 is not com- 
pletely balanced by photodeexcitation, and the resulting transi- 
tion rates are positive, that is, Y\ j + Z\j > 0. 

Populations of levels i > 3 undergo both collisional and 
radiative deexcitation. The net rates for transitions 2-3 and 
2-4 are shown in panel (d) in Fig. [6] Thus, the formation of 
the Balmer line radiation takes place at the rear of the shock 
wave in the hydrogen recombination zone in layers around the 
minimum of Z^j < 0. For the model shown in Fig.[6]the con- 
dition Ff (Ha) = corresponds to a hydrogen ionization degree 
of xu w 0.87, whereas for H/3 and Hy frequency intervals the 
condition of zero flux is fulfilled at xu ~ 0.74. 

Transfer of the Balmer line radiation in the shock wave is 
illustrated in Fig. [7] where the monochromatic radiation flux 
F v in Ha, H/3 and Hy frequency intervals is shown in three- 
dimensional plots as a function of spatial coordinate X and fre- 
quency v. The plots are given on a linear scale with respect 
to X and for the sake of convenience we discarded the outer- 
most preshock zones where the flux F v is almost independent 
of X. The frequency dependence of the flux is expressed in 
km s with negative values corresponding to the blue shifts. 
The Balmer line-center flux is about two orders of magnitude 
greater than the continuum flux, so that the dependence of the 
continuum flux on the spatial coordinate X is not seen in these 
plots. According to notations adopted the upstream flux is neg- 
ative. The frequency dependence of the space coordinate of the 
zero monochromatic flux X(F V = 0) is shown by thick contour 
lines. 



5. The emergent flux at Balmer line frequencies 

Of most interest for a comparison with the observations is the 
spectral line monochromatic radiation flux F v \ emerging from 
the upstream face of the shock wave model. In Fig. [3] we show 
the plots of F v \ as a function of frequency v for Ha, H/3 and 
Hy intervals. For the sake of convenience the flux is normal- 
ized with respect to the amplitude of the profile and is given 
with the opposite sign. It should be noted that in order to treat 
correctly the angular and frequency dependencies of the spe- 
cific intensity I v (p), the line frequency intervals considered in 
our calculations were several times of those shown in Fig. [8] 

The most conspicuous feature seen in Fig.[8]is an unshifted 
secondary maximum or a hump in the Ha and H/3 profiles. In 
the Ha frequency interval the hump is present at gas densities 
of p\ > 10" gmcrrT 3 for all values of the shock wave velocity 
considered in our study (40 kms~' < U\ < 90 kms~'), while in 
the H/3 frequency interval the weak hump appears only at p\ > 
10~ 10 gm crrT 3 and U\ < 60 km s _1 . The hump becomes less 
prominent with decreasing gas density p\ and at low densities is 
revealed as an asymmetry of the redward branch of the profile. 



log Pl = -9 T,=3000K U, = 60 km s~ 



2 - 



log /), = -!() T! = 3000K U, = 60 km s 




log p,= -12 Ti = 3000K U, = 60 km s 




-50 
km s" 



Fig. 9. The monochromatic flux emerging from the upstream 
face of the slab F v \ in the Ha frequency interval for the moving 
(solid lines) and for the static (dashed lines) medium for shock 
wave models with 10~ 12 gm crrT 3 < p\ < 10~ 9 gm crrT 3 and 
Ui = 6010ns- 1 . 



The hump in the Ha and H/3 flux profiles results from the 
radiation emitted by the narrow layer of the gas just ahead of 
the discontinuous jump. The contribution from the preshock 
gas is clearly seen when comparing the emergent flux profiles 
calculated for the moving and for the static medium. As is seen 
in Figs.|9]to[n the frequency dependence of the emergent Ha 
flux for the static medium with p\ > 10~ 1() gmcrrT 3 can be rep- 
resented as the superposition of the narrow and the broad emis- 
sion profiles with origin in the preshock and postshock regions, 
respectively. The preshock emission nearly disappears for the 
line-center optical thickness of the shock wave of r < 1. 

For U\ ^ 70 kms -1 the maximum of the preshock emission 
is blueshifted by one or a few kms -1 . This is due to the fact that 
the growth of the gas temperature in the gas element approach- 
ing the discontinuous jump is accompanied by an increase of 
the gas density p and, therefore, by a decrease of the gas flow 
velocity U. Thus, the preshock gas ahead of the discontinu- 
ous jump is slightly accelerated toward the observer and this 
effect is most conspicuous at higher densities due to stronger 
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Fig. 10. The same as in Fig.[9]ibr shock wave models with U\ = 
70kms _1 . 



absorption of the Lyman continuum radiation. The Doppler 
shift of the redward emission component increases with in- 
creasing shock wave velocity U\ (see plots for logpi = -9 
in Figs.l9ltoll 1> and for shock wave velocities within the range 
70 km s _1 < U\ < 90 km s _1 the gas flow velocity in the 
frame of the observer just ahead of the discontinuous jump is 
1 kms" 1 < -V" < 3kms^. 

As was shown above, the preshock Balmer line emission 
zone is confined to the narrow layer ahead of the discontinu- 
ous jump where photodeexcitation rates rapidly increase while 
the gas element approaches the discontinuous jump (see the 
lower panel in Fig.0}. The hump in the flux profile disappears 
when the contribution of preshock radiative deexcitation be- 
comes small in comparison with the line flux. For Ha line the 
condition of disappearance of the redward emission feature cor- 
responds to 5F(Ha)/,F ((Ha) < 0.1, where 5F(Ha) is given by 
Eq. Jl 3i . The fraction of the preshock radiation in the line pro- 
file diminishes not only with decreasing gas density p\ but also 
with increasing shock wave velocity U\. In the latter case this 
is due to the increasing width of the profile accompanying the 
growth of the postshock temperature of heavy particles T%. 
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Fig. 11. The same as in Fig.[9]for shock wave models with U\ = 
SOkms- 1 . 



To an accuracy of a few per cent the emergent line flux 
F[(Ha) and the FWHM of the line profile are the same for the 
moving and for the static medium, respectively. Because of the 
presence of the unshifted emission feature in the profiles of the 
emergent Ha radiation flux any estimate of the Doppler shift 
6V becomes ambiguous. In the present paper we evaluated 6V 
using the same procedure for all profiles both with unshifted 
emission component and without it. To this end we determined 
the Doppler shift of the profile 5V as the frequency interval 
between the profile center at its half height and the line-center 
frequency vn. Thus, values of 5V might be underestimated at 
gas densities of p\ > 10~ 10 gmcrrT 3 because of the presence of 
the unshifted redward component. 

The parameters of the shock wave models are summarised 
in Table^where we give Doppler shifts 6V, FWHM line widths 
w and optical depths at line-center frequency t for Ha, Hfi and 
Hy frequency intervals. Gas flow velocities as well as Doppler 
shifts and line widths are expressed in units of km s _1 . In the 
first three columns of Table [2 are listed the unperturbed gas 
density p\ , the shock wave velocity U \ and the upstream Mach 
number M\. The difference of the gas flow velocities in the 
frame of the observer just ahead and just behind the discontin- 
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uous jump V~ - V + gives the amplitude of the discontinuity. 
The columns labelled V* and V M give the velocity in the frame 
of the observer for the Ha emitting layer and the limiting gas 
flow velocity corresponding to isothermal compression of the 
shocked gas. The two last columns list the values of the com- 
pression ratio at the Ha emitting layer rj* = p* lp\ and of the 
limiting isothermal compression ratio 



T,=3000K U, = 70 km s" 



r]cc = poo I pi = yM 1 



(15) 



where y = 5/3 is the adiabatic exponent of the monoatomic 
gas. 

Within the considered range of shock wave velocities 
40 km s _1 < C/i < 90 (i.e. for Mach numbers ranging within 
6.2 < Mi < 14) both Doppler shifts 6V and FWHM line widths 
w increase almost linearly with increasing shock wave velocity 
U\ . At the same time 5V and w decrease with decreasing den- 
sity of the unperturbed gas p\ . 

As is seen in Table ^ the Doppler shifts are roughly one 
third of the shock wave velocity: 6V ~ lUi. There are two 
reasons for this. First, most of the Balmer line radiation is 
emitted at the rear of the shock wave where the gas under- 
goes strong compression and, therefore, slows down appre- 
ciably. In the upper panel in Fig. ^] we show the plots of 
the postshock compression ratio pjp\ for shock wave mod- 
els with U\ = 70 km s . Filled circles in the curves indi- 
cate layers with F({Ha) = where the compression ratio is 
rj* = p*/pi ~ 20 for 10~ 10 gmcirT 3 < pi < 10~ 9 gmcirT 3 and 
rf m 63 for p\ = 10~ n gm cirT 3 . Substituting these values of 
rf into Eq. we find the gas flow velocity in the layers where 
Balmer line radiation is formed to be: -V* = 36 km s for 

gm crrT 

Thus, in the Ha emitting layer the gas 



10~ 10 gm cirT 3 < pi < 10~ 9 gm cm" 3 and -V* = 34 km s" 1 



forpi 



10 11 gmcm 3 



flow velocity in the frame of the observer is roughly one half of 
the shock wave velocity. Large values of the compression ratio 
rf - p*/pi for pi < 10" 10 gm cirT 3 are due to the small opti- 
cal thickness at the Ha central frequency and, therefore, more 
rapid radiative cooling of the gas. 

A consequence of the rapid postshock radiative cooling 
is that the gas flow velocity in the Ha- emitting layer almost 
reaches its limiting value. Indeed, the upper limit of the fi- 
nal postshock compression ratio is given by the approxima- 
tion of the isothermal shock dl5> . For low ionization degrees 
(jch <k 1) the sound speed is almost independent of the gas 
density, so that for T\ = 3000 K and U\ = 70 km s _1 we have 
P<x>l pi = 197. Substituting this value of the compression ratio 
into Eq. we obtain a lower limit of the postshock velocity 
in the frame of the observer of -V«, = 33 km s _1 , whereas the 
gas flow velocity is -V* ~ 34 km s . Plots of the postshock 
gas flow velocity in the frame of the observer are shown in the 
lower panel in Fig.^] However, the ratio of the Doppler shift 
to the gas flow velocity is 5V/V* « 0.7, i.e. the low gas flow 
velocity in the hydrogen recombination zone does not account 
for the small values of the Doppler shifts. 

The second reason for the difference between 6V and V* is 
that the radiation field produced by the shock wave is character- 
ized by an anisotropy typical for slab geometry. In Fig.ll3lwe 
show the polar diagrams of the specific intensity I v (ju) for three 
points of the Ha frequency interval with displacements from 
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Fig. 12. The postshock compression ratio (upper panel) and the 
postshock gas flow velocity in the frame of the observer (lower 
panel) in shock waves with Ui = 70 km s _1 . Solid lines: pi = 
10" 9 gmcm -3 ; dotted lines: pi = 10~ 10 gmcm 3 , dashed lines: 
pi = 10~ n gm cirT 3 . Filled circles indicate the layers of the 
formation of Balmer line radiation. 



the line center of SV = 0, 20 and 30 km s _1 . For the sake of 
simplicity the polar diagrams represent the solution of the static 
transfer equation, and the specific intensity I v (p) is expressed 
in units of the mean intensity J v . The radius vector connecting 
the origin with the point of the curve is proportional to I v {p), 
whereas the angle between the radius vector and the horizontal 
axis is the directional angle 9 = cos~ l {p). Polar diagrams in the 
upper and lower panels of Fig. [O]correspond to layers of the 
hydrogen recombination zone with Fr < (jch = 0.99) and 
Fr > (xn = 0.27), respectively. For the model represented in 
Fig.lHlthe condition F(Ha) = is fulfilled at x n = 0.87. 

For the isotropic radiation field the polar diagram is the cir- 
cle of unit radius shown in Fig.[0]by dashed lines. Deviations 
from an isotropic radiation field are quite small near the line 
center (6V ^ 10 km s -1 ) but rapidly increase for 6V ^ 
20 km s~' . Thus, in calculating the monochromatic flux 



(p)fidfi 



(16) 
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Table 1. Doppler shifts 6V, FWHM line widths w and line-center optical depths t for Ha, Wj3 and Hy frequency intervals. The 
columns labelled V~, V + , V* give the gas flow velocity in the frame of the observer just ahead of the discontinuous jump, just 
behind the discontinuous jump and at the Ha emitting layer. and rj^ are the isothermal limits of the gas flow velocity in the 
frame of the observer and the compression ratio. All models were computed for T\ = 3000 K. 
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we have a larger contribution of rays with fi + ± 1 for 6V ^ 
20 km s _1 characterized by smaller values of the projection of 
the velocity vector onto the normal. 

6. Discussion 

In our attempts to solve the transfer equation for the shock wave 
structure in the co-moving frame we encountered a severe dif- 
ficulty because of the numerical instability arising at the veloc- 
ity discontinuity. In the present study this difficulty could be 
circumvented because the role of the radiation field in spectral 
lines is quite small in comparison with that of the continuum. 
According to our estimates the Doppler shifts in the hydrogen 
lines do not affect perceptibly either the structure or the radia- 
tive losses of the shock wave. In particular, the emergent flux 
integrated over the line frequency interval was found to be the 
same within < 1 % for both the static and the moving medium. 
This allowed us to leave out the effects of the Doppler shifts 
from the global iteration procedure and to solve the transfer 
equation in the frame of the observer only in the final iteration. 

The most remarkable result of our study is that the shock 
wave models show the double emission structure in the Ha and 
HjS profiles of the emergent radiation flux which is well known 
from high resolution high signal-to-noise ratio spectroscopy. 
We showed that the redward emission feature results from the 
narrow layer just ahead of the discontinuous jump within which 
the growth of collisional ionization is accompanied by pho- 
todeexcitation onto the second atomic level. The contribution 
of preshock photodeexcitation decreases with decreasing den- 
sity of the ambient gas, therefore the double emission feature 
can be considered as a tool for diagnostics of stellar atmo- 
spheres with propagating shock waves. Here one should bear 
in mind, however, that in the framework of our model the un- 
pertubed gas is at rest with respect to the observer, whereas 
in pulsating stars the gas ahead of the shock wave falls down 



onto the star with velocity in the range of one to a few dozen 
km s _1 . Thus, the preshock emission feature in Ha and H{3 pro- 
files should be observed as a redshifted component. 

Another important conclusion is that the velocity inferred 
from Doppler shifts of Balmer lines is roughly one-third of the 
shock wave velocity: 6V « | U\ . This is due to the fact that 
the gas layers emitting the Balmer line radiation are located 
at the rear of the shock wave in the hydrogen recombination 
zone where the velocity in the frame of the observer is roughly 
one half of the shock wave velocity: -V* » \U\. The ratio of 
the Doppler shift to the gas flow velocity of 6V/V* « 0.7 re- 
sults from the small optical thickness of the shock wave model 
and the anisotropy of the radiation field produced by the shock 
wave. 

Ha is the broadest emission line with FWHM comparable 
to the velocity of the shock wave U i . This, as well as the sig- 
nificant contribution from the preshock photodeexcitation zone 
make this emission line less appropriate for inferring the veloc- 
ity of the shock wave from observationally measured Doppler 
shifts. However the width of the emission profile is propor- 
tional to the temperature of the hydrogen atoms behind the dis- 
continuous jump which is related to the shock wave velocity 
U\ via the Rankine-Hugoniot relations. Thus, the width of the 
Balmer line is a function of two general quantities: the shock 
wave velocity U\ and the ambient gas density p\. In particular, 
the FWHM of Ha decreases by a factor of two with decreasing 
gas density within 10~ 9 gmcrrT 3 < p\ < 10~ 12 gmcnT 3 . 

Though our results are consistent with the observations, 
there is a number of other parameters that determine the propa- 
gation of the shock wave in the atmospheres of pulsating stars. 
In particular, our model is confined to a flat finite slab and 
therefore does not take into account the cool hydrogen gas of 
the outer stellar atmosphere above the propagating shock wave. 
Thus, we cannot exclude the role of absorption in the formation 
of the double emission structure in Balmer lines. The presence 



Yu.A. Fadeyev & D. Gillet: The structure of radiative shock waves. V 



11 



log Pl =-io T 1= 3000K U 1= 70 km s -» Bidelman, W. P. & Ratcliffe, A. E. 1954, PASP, 66, 255 

Callaway, J. 1994, ADNDT, 57, 9 

Fadeyev, Y. A. & Gillet, D. 1998, A&A, 333, 687 jPaner J) 
Fadeyev, Y. A. & Gillet, D. 2000, A&A, 354, 349 jPaner iA 
Fadeyev, Y. A. & Gillet, D. 2001, A&A, 368, 901 JPaper IlJ) 
Fadeye v, Y. A., L e Coroller, H., & Gillet, D. 2002, A&A, 392, 
735 dPaner IVh 

Fox, M. W., Wood, P. R., & Dopita, M. A. 1984, ApJ, 286, 337 
Gillet, D. 1988, A&A, 192, 206 
Gorbatskii, V. G. 1961, Soviet Astronomy, 5, 192 
Hutchinson, J. L., Lillie, C. F, & Hill, S. J. 1975, AJ, 80, 1044 
Joy, A. H. 1947, ApJ, 106, 288 
Joy, A. H. 1954, ApJS, 1,39 

Kraft, R. P., Camp, D. C, & Hughes, W. T. 1959, ApJ, 130, 90 
Lebre, A. & Gillet, D. 1991, A&A, 246, 490 
Lebre, A. & Gillet, D. 1992, A&A, 255, 221 
Lennon, D. J., Brown, P. J. F, Dufton, P. L., & Lynas-Gray, 

A. E. 1986, MNRAS, 222, 719 
Mihalas, D. 1978, Stellar atmospheres /2nd edition/ (San 

Francisco, W. H. Freeman and Co., 1978. 650 p.) 
Mihalas, D. 1980, ApJ, 238, 1042 

Mihalas, D., Kunasz, P. B., & Hummer, D. G. 1975, ApJ, 202, 
465 

Preston, G. W. 1962, ApJ, 136, 866 
Preston, G. W. 1964, ARA&A, 2, 23 
Richter, H. & Wood, P. R. 2001, A&A, 369, 1027 
Rybicki, G. B. & Hummer, D. G. 1991, A&A, 245, 171 
Scholz, T. T., Walters, H. R. J., Burke, P. J., & Scott, M. P. 1990, 

MNRAS, 242, 692 
Wallerstein, G. 1959, ApJ, 130, 560 
Willson, L. A. 1976, ApJ, 205, 172 
Woodsworth, A. W. 1995, ApJ, 444, 396 
-1 o 1 

Fig. 13. Polar diagrams of the normalized specific intensity 
I v (fi) in the frequency interval of the Ha line in the hydrogen 
recombination zone. Upper panel: layers with F v < 0, lower 
panel: layers with F v > 0. Plots of I v (/j)/J v represent the static 
solution of the transfer equation. The values of AV give the 
distance from the line-center frequency vo expressed in units 
of km s . Dashed line: the polar diagram of the isotropic radi- 
ation field. 

of such absorption in emission profiles observed in R V Tau and 
W Vir stars was pointed out, for example, by iLebre & Gillet! 
ill99lHl992h . 
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Fig. 7. The monochromatic flux F Y in the frequency intervals of the Ha, Hfi and Hy lines. Thick lines show the level of F v = 0. 




Fig. 8. The normalized monochromatic flux at the upstream face of the slab F v \ in the Ha, R/3 and Ry frequency intervals for 
shock wave models withpi = 10" 11 , lO" 10 , 10~ 9 gm cm" 3 , T\ = 3000 K, Ui = 50, 70, 90 km s" 1 . 



